Priming With Rhinovirus Protects Mice Against a Lethal Pulmonary Coronavirus Infection.

Rhinoviruses (RV) have been shown to inhibit subsequent infection by heterologous respiratory viruses, including influenza viruses and severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2). To better understand the mechanisms whereby RV protects against pulmonary coronavirus infection, we used a native murine virus, mouse hepatitis virus strain 1 (MHV-1), that causes severe disease in the lungs of infected mice. We found that priming of the respiratory tract with RV completely prevented mortality and reduced morbidity of a lethal MHV-1 infection. Replication of MHV-1 was reduced in RV-primed mouse lungs although expression of antiviral type I interferon, IFN-ß, was more robust in mice infected with MHV-1 alone. We further showed that signaling through the type I interferon receptor was required for survival of mice given a non-lethal dose of MHV-1. RV-primed mice had reduced pulmonary inflammation and hemorrhage and influx of leukocytes, especially neutrophils, in the airways upon MHV-1 infection. Although MHV-1 replication was reduced in RV-primed mice, RV did not inhibit MHV-1 replication in coinfected lung epithelial cells in vitro. In summary, RV-mediated priming in the respiratory tract reduces viral replication, inflammation, and tissue damage, and prevents mortality of a pulmonary coronavirus infection in mice. These results contribute to our understanding of how distinct respiratory viruses interact with the host to affect disease pathogenesis, which is a critical step in understanding how respiratory viral coinfections impact human health.

Global Composition of the Bacteriophage Community in Honey Bees.

The microbial communities in animal digestive systems are critical for host development and health. They stimulate the immune system during development, synthesize important chemical compounds like hormones, aid in digestion, competitively exclude pathogens, etc. Compared to the bacterial and fungal components of the microbiome, we know little about the temporal and spatial dynamics of bacteriophage communities in animal digestive systems. Recently, the bacteriophages of the honey bee gut were characterized in two European bee populations. Most of the bacteriophages described in these two reports were novel, harbored many metabolic genes in their genomes, and had a community structure that suggests coevolution with their bacterial hosts. To describe the conservation of bacteriophages in bees and begin to understand their role in the bee microbiome, we sequenced the virome of Apis mellifera from Austin, TX, and compared bacteriophage compositions among three locations around the world. We found that most bacteriophages from Austin are novel, sharing no sequence similarity with anything in public repositories. However, many bacteriophages are shared among the three bee viromes, indicating specialization of bacteriophages in the bee gut. Our study, along with the two previous bee virome studies, shows that the bee gut bacteriophage community is simple compared to that of many animals, consisting of several hundred types of bacteriophages that primarily infect four of the dominant bacterial phylotypes in the bee gut. IMPORTANCE Viruses that infect bacteria (bacteriophages) are abundant in the microbial communities that live on and in plants and animals. However, our knowledge of the structure, dynamics, and function of these viral communities lags far behind our knowledge of their bacterial hosts. We sequenced the first bacteriophage community of honey bees from the United States and compared the U.S. honey bee bacteriophage community to those of samples from Europe. Our work is an important characterization of an economically critical insect species and shows how bacteriophage communities can contain highly conserved individuals and be highly variable in composition across a wide geographic range.

Nuku, a family of primate retrocopies derived from KU70.

The gene encoding the ubiquitous DNA repair protein, Ku70p, has undergone extensive copy number expansion during primate evolution. Gene duplications of KU70 have the hallmark of long interspersed element-1 mediated retrotransposition with evidence of target-site duplications, the poly-A tails, and the absence of introns. Evolutionary analysis of this expanded family of KU70-derived "NUKU" retrocopies reveals that these genes are both ancient and also actively being created in extant primate species. NUKU retrocopies show evidence of functional divergence away from KU70, as evinced by their altered pattern of tissue expression and possible tissue-specific translation. Molecular modeling predicted that amino acid changes in Nuku2p at the interaction interface with Ku80p would prevent the assembly of the Ku heterodimer. The lack of Nuku2p-Ku80p interaction was confirmed by yeast two-hybrid assay, which contrasts the robust interaction of Ku70p-Ku80p. While several NUKU retrocopies appear to have been degraded by mutation, NUKU2 shows evidence of positive natural selection, suggesting that this retrocopy is undergoing neofunctionalization. Although Nuku proteins do not appear to antagonize retrovirus transduction in cell culture, the observed expansion and rapid evolution of NUKUs could be being driven by alternative selective pressures related to infectious disease or an undefined role in primate physiology.

Rhinovirus Reduces the Severity of Subsequent Respiratory Viral Infections by Interferon-Dependent and -Independent Mechanisms.

Coinfection by heterologous viruses in the respiratory tract is common and can alter disease severity compared to infection by individual virus strains. We previously found that inoculation of mice with rhinovirus (RV) 2 days before inoculation with a lethal dose of influenza A virus [A/Puerto Rico/8/34 (H1N1) (PR8)] provides complete protection against mortality. Here, we extended that finding to a second lethal respiratory virus, pneumonia virus of mice (PVM), and analyzed potential mechanisms of RV-induced protection. RV completely prevented mortality and weight loss associated with PVM infection. Major changes in host gene expression upon PVM infection were delayed compared to PR8. RV induced earlier recruitment of inflammatory cells, which were reduced at later times in RV-inoculated mice. Findings common to both virus pairs included the upregulated expression of mucin-associated genes and dampening of inflammation-related genes in mice that were inoculated with RV before lethal virus infection. However, type I interferon (IFN) signaling was required for RV-mediated protection against PR8 but not PVM. IFN signaling had minor effects on PR8 replication and contributed to controlling neutrophilic inflammation and hemorrhagic lung pathology in RV/PR8-infected mice. These findings, combined with differences in virus replication levels and disease severity, suggest that the suppression of inflammation in RV/PVM-infected mice may be due to early, IFN-independent suppression of viral replication, while that in RV/PR8-infected mice may be due to IFN-dependent modulation of immune responses. Thus, a mild upper respiratory viral infection can reduce the severity of a subsequent severe viral infection in the lungs through virus-dependent mechanisms. IMPORTANCE Respiratory viruses from diverse families cocirculate in human populations and are frequently detected within the same host. Although clinical studies suggest that infection by multiple different respiratory viruses may alter disease severity, animal models in which we can control the doses, timing, and strains of coinfecting viruses are critical to understanding how coinfection affects disease severity. Here, we compared gene expression and immune cell recruitment between two pairs of viruses (RV/PR8 and RV/PVM) inoculated sequentially in mice, both of which result in reduced severity compared to lethal infection by PR8 or PVM alone. Reduced disease severity was associated with suppression of inflammatory responses in the lungs. However, differences in disease kinetics and host and viral gene expression suggest that protection by coinfection with RV may be due to distinct molecular mechanisms. Indeed, we found that antiviral cytokine signaling was required for RV-mediated protection against lethal infection by PR8 but not PVM.

ΦX174 Attenuation by Whole-Genome Codon Deoptimization.

Natural selection acting on synonymous mutations in protein-coding genes influences genome composition and evolution. In viruses, introducing synonymous mutations in genes encoding structural proteins can drastically reduce viral growth, providing a means to generate potent, live-attenuated vaccine candidates. However, an improved understanding of what compositional features are under selection and how combinations of synonymous mutations affect viral growth is needed to predictably attenuate viruses and make them resistant to reversion. We systematically recoded all nonoverlapping genes of the bacteriophage ΦX174 with codons rarely used in its Escherichia coli host. The fitness of recombinant viruses decreases as additional deoptimizing mutations are made to the genome, although not always linearly, and not consistently across genes. Combining deoptimizing mutations may reduce viral fitness more or less than expected from the effect size of the constituent mutations and we point out difficulties in untangling correlated compositional features. We test our model by optimizing the same genes and find that the relationship between codon usage and fitness does not hold for optimization, suggesting that wild-type ΦX174 is at a fitness optimum. This work highlights the need to better understand how selection acts on patterns of synonymous codon usage across the genome and provides a convenient system to investigate the genetic determinants of virulence.

Plasmodium falciparum DHFR and DHPS Mutations Are Associated With HIV-1 Co-Infection and a Novel DHPS Mutation I504T Is Identified in Western Kenya.

Antifolate resistance is significant in Kenya and presumed to result from extensive use and cross-resistance between antifolate antimalarials and antibiotics, including cotrimoxazole/Bactrim used for HIV-1 chemotherapy. However, little is known about antifolate-resistant malaria in the context of newly diagnosed HIV-1 co-infection prior to administration of HIV-1 chemotherapy. Blood samples from a cross-sectional study of asymptomatic adult Kenyans enrolled during voluntary HIV testing were analyzed by PCR for Plasmodium spp. More than 95% of volunteers with identifiable parasite species (132 HIV-1 co-infected) were infected with Plasmodium falciparum alone or P. falciparum with Plasmodium ovale and/or Plasmodium malariae. Deep sequencing was used to screen for mutations in P. falciparum dihydrofolate reductase (dhfr) (N51I, C59R, S108N, I164L) and dihydropteroate synthase (dhps) (S436H, A437G, K540E, A581G) from 1133 volunteers. Individual mutations in DHPS but not DHFR correlated with HIV-1 status. DHFR haplotype diversity was significantly different among volunteers by gender and HIV-1 status. DHPS haplotype diversity by HIV-1 status was significantly different between volunteers paired by age and gender, indicating that patterns of resistance were independent of these variables. Molecular simulations for a novel DHPS mutation (I504T) suggested that the mutated protein has increased affinity for the endogenous ligand DHPPP and decreased affinity for drug binding. A sub-group of monoclonal infections revealed that age and parasitemia were not correlated and enabled identification of a rare septuple-mutant haplotype (IRNL-HGEA). In our study, adult Kenyans newly diagnosed with HIV-1 infection were predominantly infected with moderately resistant P. falciparum, with patterns of infecting parasite genotypes significantly associated with HIV-1 status. Together with the discovery of DHPS I504T, these data indicate that antifolate resistance continues to evolve in Kenya. Further, they highlight the need to understand the effects of associated mutations on both fitness and resistance of P. falciparum in the context of HIV-1 co-infection to better inform treatment for asymptomatic malaria.

Organic Farming Sharpens Plant Defenses in the Field.

Plants deploy a variety of chemical and physical defenses to protect themselves against herbivores and pathogens. Organic farming seeks to enhance these responses by improving soil quality, ultimately altering bottom up regulation of plant defenses. While laboratory studies suggest this approach is effective, it remains unclear whether organic agriculture encourages more-active plant defenses under real-world conditions. Working on the farms of cooperating growers, we examined gene expression in the leaves of two potato (Solanum tuberosum) varieties, grown on organic vs. conventional farms. For one variety, Norkotah, we found significantly heightened initiation of genes associated with plant-defense pathways in plants grown in organic vs. conventional fields. Organic Norkotah fields exhibited lower levels of nitrate in soil and of nitrogen in plant foliage, alongside differences in communities of soil bacteria, suggesting possible links between soil management and observed differences in plant defenses. Additionally, numbers of predatory and phloem-feeding insects were higher in organic than conventional fields. A second potato variety, Alturas, which is generally grown using fewer inputs and in poorer-quality soils, exhibited lower overall herbivore and predator numbers, few differences in soil ecology, and no differences in gene-activity in organic and conventional farming systems. Altogether, our results suggest that organic farming has the potential to increase plants' resistance to herbivores, possibly facilitating reduced need for insecticide applications. These benefits appear to be mediated by plant variety and/or farming context.

Saturation Mutagenesis Genome Engineering of Infective ΦX174 Bacteriophage via Unamplified Oligo Pools and Golden Gate Assembly.

Here we present a novel protocol for the construction of saturation single-site-and massive multisite-mutant libraries of a bacteriophage. We segmented the ΦX174 genome into 14 nontoxic and nonreplicative fragments compatible with Golden Gate assembly. We next used nicking mutagenesis with oligonucleotides prepared from unamplified oligo pools with individual segments as templates to prepare near-comprehensive single-site mutagenesis libraries of genes encoding the F capsid protein (421 amino acids scanned) and G spike protein (172 amino acids scanned). Libraries possessed greater than 99% of all 11 860 programmed mutations. Golden Gate cloning was then used to assemble the complete ΦX174 mutant genome and generate libraries of infective viruses. This protocol will enable reverse genetics experiments for studying viral evolution and, with some modifications, can be applied for engineering therapeutically relevant bacteriophages with larger genomes.

Organic Soils Control Beetle Survival While Competitors Limit Aphid Population Growth.

Soil chemistry and microbial diversity can impact the vigor and nutritive qualities of plants, as well as plants' ability to deploy anti-herbivore defenses. Soil qualities often vary dramatically on organic versus conventional farms, reflecting the many differences in soil management practices between these farming systems. We examined soil-mediated effects on herbivore performance by growing potato plants (Solanum tuberosum L.) in soils collected from organic or conventional commercial farm fields, and then exposing these plants to herbivory by green peach aphids (Myzus persicae Sulzer, Hemiptera: Aphididae) and/or Colorado potato beetles (Leptinotarsa decemlineata Say, Coleoptera: Chrysomelidae). Responses of the two potato pests varied dramatically. Survivorship of Colorado potato beetles was almost 3× higher on plants grown in organic than in conventional soils, but was unaffected by the presence of aphids. In contrast, aphid colony growth was twice as rapid when aphids were reared alone rather than with Colorado potato beetles, but was unaffected by soil type. We saw no obvious differences in soil nutrients when comparing organic and conventional soils. However, we saw a higher diversity of bacteria in organic soils, and potato plants grown in this soil had a lower carbon concentration in foliar tissue. In summary, the herbivore species differed in their susceptibility to soil- versus competitor-mediated effects, and these differences may be driven by microbe-mediated changes in host plant quality. Our results suggest that soil-mediated effects on pest growth can depend on herbivore species and community composition, and that soil management strategies that promote plant health may also increase host quality for pests.

Cicada Endosymbionts Have tRNAs That Are Correctly Processed Despite Having Genomes That Do Not Encode All of the tRNA Processing Machinery.

Gene loss and genome reduction are defining characteristics of endosymbiotic bacteria. The most highly reduced endosymbiont genomes have lost numerous essential genes related to core cellular processes such as replication, transcription, and translation. Computational gene predictions performed for the genomes of the two bacterial symbionts of the cicada Diceroprocta semicincta, "Candidatus Hodgkinia cicadicola" (Alphaproteobacteria) and "Ca Sulcia muelleri" (Bacteroidetes), have found only 26 and 16 tRNA genes and 15 and 10 aminoacyl tRNA synthetase genes, respectively. Furthermore, the original "Ca Hodgkinia cicadicola" genome annotation was missing several essential genes involved in tRNA processing, such as those encoding RNase P and CCA tRNA nucleotidyltransferase as well as several RNA editing enzymes required for tRNA maturation. How these cicada endosymbionts perform basic translation-related processes remains unknown. Here, by sequencing eukaryotic mRNAs and total small RNAs, we show that the limited tRNA set predicted by computational annotation of "Ca Sulcia muelleri" and "Ca Hodgkinia cicadicola" is likely correct. Furthermore, we show that despite the absence of genes encoding tRNA processing activities in the symbiont genomes, symbiont tRNAs have correctly processed 5' and 3' ends and seem to undergo nucleotide modification. Surprisingly, we found that most "Ca Hodgkinia cicadicola" and "Ca Sulcia muelleri" tRNAs exist as tRNA halves. We hypothesize that "Ca Sulcia muelleri" and "Ca Hodgkinia cicadicola" tRNAs function in bacterial translation but require host-encoded enzymes to do so.IMPORTANCE The smallest bacterial genomes, in the range of about 0.1 to 0.5 million base pairs, are commonly found in the nutritional endosymbionts of insects. These tiny genomes are missing genes that encode proteins and RNAs required for the translation of mRNAs, one of the most highly conserved and important cellular processes. In this study, we found that the bacterial endosymbionts of cicadas have genomes which encode incomplete tRNA sets and lack genes required for tRNA processing. Nevertheless, we found that endosymbiont tRNAs are correctly processed at their 5' and 3' ends and, surprisingly, that mostly exist as tRNA halves. We hypothesize that the cicada host must supply its symbionts with these missing tRNA processing activities.

Comparative genomics of Bifidobacterium species isolated from marmosets and humans.

The genus Bifidobacterium is purported to have beneficial consequences for human health and is a major component of many gastrointestinal probiotics. Although species of Bifidobacterium are generally at low relative frequency in the adult human gastrointestinal tract, they can constitute high proportions of the gastrointestinal communities of adult marmosets. To identify genes that might be important for the maintenance of Bifidobacterium in adult marmosets, ten strains of Bifidobacterium were isolated from the feces of seven adult marmosets, and their genomes were sequenced. There were six B. reuteri strains, two B. callitrichos strains, one B. myosotis sp. nov. and one B. tissieri sp. nov. among our isolates. Phylogenetic analysis showed that three of the four species we isolated were most closely related to B. bifidum, B. breve and B. longum, which are species found in high abundance in human infants. There were 1357 genes that were shared by at least one strain of B. reuteri, B. callitrichos, B. breve, and B. longum, and 987 genes that were found in all strains of the four species. There were 106 genes found in B. reuteri and B. callitrichos but not in human bifidobacteria, and several of these genes were involved in nutrient uptake. These pathways for nutrient uptake appeared to be specific to Bifidobacterium from New World monkeys. Additionally, the distribution of Bifidobacterium in fecal samples from captive adult marmosets constituted as much as 80% of the gut microbiome, although this was variable between individuals and colonies. We suggest that nutrient transporters may be important for the maintenance of Bifidobacterium during adulthood in marmosets.

A Rapid Method for Sequencing Double-Stranded RNAs Purified from Yeasts and the Identification of a Potent K1 Killer Toxin Isolated from Saccharomyces cerevisiae.

Mycoviruses infect a large number of diverse fungal species, but considering their prevalence, relatively few high-quality genome sequences have been determined. Many mycoviruses have linear double-stranded RNA genomes, which makes it technically challenging to ascertain their nucleotide sequence using conventional sequencing methods. Different specialist methodologies have been developed for the extraction of double-stranded RNAs from fungi and the subsequent synthesis of cDNAs for cloning and sequencing. However, these methods are often labor-intensive, time-consuming, and can require several days to produce cDNAs from double-stranded RNAs. Here, we describe a comprehensive method for the rapid extraction and sequencing of dsRNAs derived from yeasts, using short-read next generation sequencing. This method optimizes the extraction of high-quality double-stranded RNAs from yeasts and 3' polyadenylation for the initiation of cDNA synthesis for next-generation sequencing. We have used this method to determine the sequence of two mycoviruses and a double-stranded RNA satellite present within a single strain of the model yeast Saccharomyces cerevisiae. The quality and depth of coverage was sufficient to detect fixed and polymorphic mutations within viral populations extracted from a clonal yeast population. This method was also able to identify two fixed mutations within the alpha-domain of a variant K1 killer toxin encoded on a satellite double-stranded RNA. Relative to the canonical K1 toxin, these newly reported mutations increased the cytotoxicity of the K1 toxin against a specific species of yeast.

One Hundred Mitochondrial Genomes of Cicadas.

Mitochondrial genomes can provide valuable information on the biology and evolutionary histories of their host organisms. Here, we present and characterize the complete coding regions of 107 mitochondrial genomes (mitogenomes) of cicadas (Insecta: Hemiptera: Auchenorrhyncha: Cicadoidea), representing 31 genera, 61 species, and 83 populations. We show that all cicada mitogenomes retain the organization and gene contents thought to be ancestral in insects, with some variability among cicada clades in the length of a region between the genes nad2 and cox1, which encodes 3 tRNAs. Phylogenetic analyses using these mitogenomes recapitulate a recent 5-gene classification of cicadas into families and subfamilies, but also identify a species that falls outside of the established taxonomic framework. While protein-coding genes are under strong purifying selection, tests of relative evolutionary rates reveal significant variation in evolutionary rates across taxa, highlighting the dynamic nature of mitochondrial genome evolution in cicadas. These data will serve as a useful reference for future research into the systematics, ecology, and evolution of the superfamily Cicadoidea.

Multiple origins of interdependent endosymbiotic complexes in a genus of cicadas.

Bacterial endosymbionts that provide nutrients to hosts often have genomes that are extremely stable in structure and gene content. In contrast, the genome of the endosymbiont Hodgkinia cicadicola has fractured into multiple distinct lineages in some species of the cicada genus Tettigades To better understand the frequency, timing, and outcomes of Hodgkinia lineage splitting throughout this cicada genus, we sampled cicadas over three field seasons in Chile and performed genomics and microscopy on representative samples. We found that a single ancestral Hodgkinia lineage has split at least six independent times in Tettigades over the last 4 million years, resulting in complexes of between two and six distinct Hodgkinia lineages per host. Individual genomes in these symbiotic complexes differ dramatically in relative abundance, genome size, organization, and gene content. Each Hodgkinia lineage retains a small set of core genes involved in genetic information processing, but the high level of gene loss experienced by all genomes suggests that extensive sharing of gene products among symbiont cells must occur. In total, Hodgkinia complexes that consist of multiple lineages encode nearly complete sets of genes present on the ancestral single lineage and presumably perform the same functions as symbionts that have not undergone splitting. However, differences in the timing of the splits, along with dissimilar gene loss patterns on the resulting genomes, have led to very different outcomes of lineage splitting in extant cicadas.

Selecting among three basic fitness landscape models: Additive, multiplicative and stickbreaking.

Fitness landscapes map genotypes to organismal fitness. Their topographies depend on how mutational effects interact - epistasis - andare important for understanding evolutionary processes such as speciation, the rate of adaptation, the advantage of recombination, and the predictability versus stochasticity of evolution. The growing amount of data has made it possible to better test landscape models empirically. We argue that this endeavor will benefit from the development and use of meaningful basic models against which to compare more complex models. Here we develop statistical and computational methods for fitting fitness data from mutation combinatorial networks to three simple models: additive, multiplicative and stickbreaking. We employ a Bayesian framework for doing model selection. Using simulations, we demonstrate that our methods work and we explore their statistical performance: bias, error, and the power to discriminate among models. We then illustrate our approach and its flexibility by analyzing several previously published datasets. An R-package that implements our methods is available in the CRAN repository under the name Stickbreaker.

Modeling time-series data from microbial communities.

As sequencing technologies have advanced, the amount of information regarding the composition of bacterial communities from various environments (for example, skin or soil) has grown exponentially. To date, most work has focused on cataloging taxa present in samples and determining whether the distribution of taxa shifts with exogenous covariates. However, important questions regarding how taxa interact with each other and their environment remain open thus preventing in-depth ecological understanding of microbiomes. Time-series data from 16S rDNA amplicon sequencing are becoming more common within microbial ecology, but methods to infer ecological interactions from these longitudinal data are limited. We address this gap by presenting a method of analysis using Poisson regression fit with an elastic-net penalty that (1) takes advantage of the fact that the data are time series; (2) constrains estimates to allow for the possibility of many more interactions than data; and (3) is scalable enough to handle data consisting of thousands of taxa. We test the method on gut microbiome data from white-throated woodrats (Neotoma albigula) that were fed varying amounts of the plant secondary compound oxalate over a period of 22 days to estimate interactions between OTUs and their environment.

Genome expansion via lineage splitting and genome reduction in the cicada endosymbiont Hodgkinia.

Comparative genomics from mitochondria, plastids, and mutualistic endosymbiotic bacteria has shown that the stable establishment of a bacterium in a host cell results in genome reduction. Although many highly reduced genomes from endosymbiotic bacteria are stable in gene content and genome structure, organelle genomes are sometimes characterized by dramatic structural diversity. Previous results from Candidatus Hodgkinia cicadicola, an endosymbiont of cicadas, revealed that some lineages of this bacterium had split into two new cytologically distinct yet genetically interdependent species. It was hypothesized that the long life cycle of cicadas in part enabled this unusual lineage-splitting event. Here we test this hypothesis by investigating the structure of the Ca. Hodgkinia genome in one of the longest-lived cicadas, Magicicada tredecim. We show that the Ca. Hodgkinia genome from M. tredecim has fragmented into multiple new chromosomes or genomes, with at least some remaining partitioned into discrete cells. We also show that this lineage-splitting process has resulted in a complex of Ca. Hodgkinia genomes that are 1.1-Mb pairs in length when considered together, an almost 10-fold increase in size from the hypothetical single-genome ancestor. These results parallel some examples of genome fragmentation and expansion in organelles, although the mechanisms that give rise to these extreme genome instabilities are likely different.

Biochemical association of metabolic profile and microbiome in chronic pressure ulcer wounds.

Chronic, non-healing wounds contribute significantly to the suffering of patients with co-morbidities in the clinical population with mild to severely compromised immune systems. Normal wound healing proceeds through a well-described process. However, in chronic wounds this process seems to become dysregulated at the transition between resolution of inflammation and re-epithelialization. Bioburden in the form of colonizing bacteria is a major contributor to the delayed headlining in chronic wounds such as pressure ulcers. However how the microbiome influences the wound metabolic landscape is unknown. Here, we have used a Systems Biology approach to determine the biochemical associations between the taxonomic and metabolomic profiles of wounds colonized by bacteria. Pressure ulcer biopsies were harvested from primary chronic wounds and bisected into top and bottom sections prior to analysis of microbiome by pyrosequencing and analysis of metabolome using 1H nuclear magnetic resonance (NMR) spectroscopy. Bacterial taxonomy revealed that wounds were colonized predominantly by three main phyla, but differed significantly at the genus level. While taxonomic profiles demonstrated significant variability between wounds, metabolic profiles shared significant similarity based on the depth of the wound biopsy. Biochemical association between taxonomy and metabolic landscape indicated significant wound-to-wound similarity in metabolite enrichment sets and metabolic pathway impacts, especially with regard to amino acid metabolism. To our knowledge, this is the first demonstration of a statistically robust correlation between bacterial colonization and metabolic landscape within the chronic wound environment.

Sympatric speciation in a bacterial endosymbiont results in two genomes with the functionality of one.

Mutualisms that become evolutionarily stable give rise to organismal interdependencies. Some insects have developed intracellular associations with communities of bacteria, where the interdependencies are manifest in patterns of complementary gene loss and retention among members of the symbiosis. Here, using comparative genomics and microscopy, we show that a three-member symbiotic community has become a four-way assemblage through a novel bacterial lineage-splitting event. In some but not all cicada species of the genus Tettigades, the endosymbiont Candidatus Hodgkinia cicadicola has split into two new cytologically distinct but metabolically interdependent species. Although these new bacterial genomes are partitioned into discrete cell types, the intergenome patterns of gene loss and retention are almost perfectly complementary. These results defy easy classification: they show genomic patterns consistent with those observed after both speciation and whole-genome duplication. We suggest that our results highlight the potential power of nonadaptive forces in shaping organismal complexity.

Dynamic recruitment of amino acid transporters to the insect/symbiont interface.

Symbiosis is well known to influence bacterial symbiont genome evolution and has recently been shown to shape eukaryotic host genomes. Intriguing patterns of host genome evolution, including remarkable numbers of gene duplications, have been observed in the pea aphid, a sap-feeding insect that relies on a bacterial endosymbiont for amino acid provisioning. Previously, we proposed that gene duplication has been important for the evolution of symbiosis based on aphid-specific gene duplication in amino acid transporters (AATs), with some paralogs highly expressed in the cells housing symbionts (bacteriocytes). Here, we use a comparative approach to test the role of gene duplication in enabling recruitment of AATs to bacteriocytes. Using genomic and transcriptomic data, we annotate AATs from sap-feeding and non sap-feeding insects and find that, like aphids, AAT gene families have undergone independent large-scale gene duplications in three of four additional sap-feeding insects. RNA-seq differential expression data indicate that, like aphids, the sap-feeding citrus mealybug possesses several lineage-specific bacteriocyte-enriched paralogs. Further, differential expression data combined with quantitative PCR support independent evolution of bacteriocyte enrichment in sap-feeding insect AATs. Although these data indicate that gene duplication is not necessary to initiate host/symbiont amino acid exchange, they support a role for gene duplication in enabling AATs to mediate novel host/symbiont interactions broadly in the sap-feeding suborder Sternorrhyncha. In combination with recent studies on other symbiotic systems, gene duplication is emerging as a general pattern in host genome evolution.